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ABSTRACT 


With  the  use  of  moving  instruments,  such  as  airborne  gravimeters,  which 
are  often  related  to  precise  inertial  stabilization  systems,  the  essentially  statical 
character  of  geodesy  begins  to  be  enlarged  by  kinematical  features.  The  potential 
of  gravity,  combining  gravitational  attraction  and  centrifugal  force,  is  no  longer 
adequate  for  kinematics,  but  other  inertial  forces  of  Coriolis  type  must  also  be 
considered. 


The  main  purpose  of  the  present  paper  is  an  investigation  of  the  geodetic 
aspects  of  the  interrelation  of  gravitational  and  inertial  forces  and  their  separa¬ 
tion  by  means  of  structural  differences  in  their  respective  fields.  For  a  deeper 
insight,  the  general  theory  of  relativity  is  indispensable;  it  also  furnishes  con¬ 
venient  mathematical  techniques. 

The  extraction  of  purely  gravitational  effects  is  possible  with  second  and 
third  derivatives  of  the  potential;  therefore  geodetic  applications  of  these  quantities 
are  discussed,  integral  formulas  similar  to  Stokes'  integral  being  given  for  this 
purpose. 
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KINEMATICAL  GEODESY 


Introduction 

Until  recently,  geodetic  measurements  were  visually  performed  with 
stationary  instruments,  that  is,  instruments  at  rest  on  the  earth's  surface, 
such  as  theodolites  or  gravimeters.  Therefore,  only  the  statical  features  of 
the  earth  and  its  gravity  field  were  of  geodetic  interest.  A  typical  example  is 
the  omission  of  the  Coriolis  force  from  the  gravity  field  of  the  earth,  because 
this  force  acts  on  moving  bodies  only. 

With  the  advent  of  artificial  satellites  and  precise  inertial  navigation 
systems,  and  with  the  development  of  shipborne  and  airborne  systems  for 
gravity  measurement,  the  element  of  time  is  now  entering  into  physical  geodesy; 
the  essentially  statical  character  of  geodesy  begins  to  be  enlarged  by  kinematical 
features.  As  a  simple  example,  Coriolis  forces  show  up  in  marine  and  aerial 
gravimetry.  This  indicates  that  th^  potential  of  gravity,  combining  gravitational 
attraction  and  centrifugal  force  but  excluding  the  Coriolis  force,  is  no  longer 
adequate  for  treating  certain  kinematical  effects  of  the  earth's  gravitational  field. 

The  Principle  of  Equivalence,  as  expressed  by  the  identity  of  inertial  and 
gravitational  mass,  makes  it  possible  to  combine  gravitational  attraction  and 
centrifugal  force,  which  is  an  inertial  effect,  into  a  unified  force,  gravity.  The 
same  principle,  however,  also  makes  it  extremely  difficult  to  separate  the  effect 
of  undesirable  inertial  forces  from  the  genuine  gravitational  forces,  in  which  we 
are  interested  in  gravity  measurement. 


This  problem  of  separation  of  inertial  and  gravitational  forces  is  important 
if  the  measuring  system  is  exposed  to  irregular  accelerations,  such  as  aboard  a 
ship  or  an  aircraft.  It  does  not  occur  in  considering  the  motion  of  a  satellite, 
f.ince  the  main  purpose  of  the  present  report  is  an  investigation  of  the  geodetic 
ispects  of  the  interrelation  of  gravitation  and  inertia,  satellite  dynamics  is 
niturally  excluded  here. 

In  a  sense,  the  problems  of  aerial  gravimetry  and  inertial  positioning 
(navigation)  are  inverse  to  each  other.  In  aerial  gravimetry,  outside  geodetic 
information  is  used  to  separate  gravitational  from  inertial  effects  in  order  to 
measure  gravity  (or  other  quantities  c‘  the  gravitational  field).  In  inertial 
navigation,  the  gravitational  field  is  needed  as  an  input  to  obtain  position  (which, 
applied  to  geodesy,  means  geodetic  position) .  Schematically  we  may  write 


Aerial  Gravimetry 


Geodetic  Position 


Gravity  Field; 
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A.  Statistical  frequency  behaviour.  The  disturbing  accelerations  usually 
have  a  higher  frequency  than  gravity;  thus  they  may  be  removed,  to  a 
certain  extent,  by  filtering  techniques. 

B.  Structural  differences  between  the  gravitational  and  the  inertial  fields. 
The  second  field  is  more  regular  than  the  first. 

Filtering  techniques  based  on  statistical  behaviour  are  discussed  in  another 
report  (Moritz,  1967).  Here  we  shall  investigate  the  mathematical  structure  of 
the  fields  involved  and  geodetic  applications. 

It  turns  out  that  factor  B  has  no  effect  on  the  gravity  g  (or  more  generally 
speaking,  on  the  first  derivatives  of  the  potential),  so  that  factor  A  must  be  used 
with  measurements  of  g.  It  is  only  in  second  and  higher  derivatives  that  the 
structural  differences  between  gravitational  and  inertial  fields  will  contribute  to 
their  separation. 

For  this  reason,  we  shall  in  Chapter  1  give  a  theoretical  (perhaps  some¬ 
what  unconventional)  outline  of  possible  geodetic  uses  of  second  and  third  deri¬ 
vatives  of  the  potential. 

Chapter  2  is  devoted  to  theoretical  aspects  of  measurement  of  these  deriva¬ 
tives,  with  emphasis  on  separability  of  inertial  and  gravitational  effects  because  of 
the  structure  of  their  respective  fields.  The  approach  through  classical  mechan  cs 
is  simple  and  straightforward.  However,  since  the  Principle  of  Equivalence  can¬ 
not  be  explained  on  the  basis  of  classical  mechanics,  but  can  be  understood  only 
on  the  basis  of  Einstein's  general  theory  of  relativity,  there  will  always  be  the 
shadow  of  a  doubt  in  arguments  involving  classical  mechanics  only;  we  shall  try 
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to  remove  this  shadow  by  a  rather  detailed  study  of  the  relativistic  aspects  of  the 


problem. 

Time  is  obviously  the  key  concept  in  the  kinematical1  applications  of  geodesy. 
Since  time  is  a  natural  part  of  the  space-time  of  the  theory  of  relativity,  we  may 
expect  conceptual  and  mathematical  advantages  in  applying  this  theory  to  kine¬ 
matical  geodesy.  As  an  example,  certain  relativistic  quantities  generalize  the 
potential  so  as  to  contain  the  Coriolis  force  as  well  as  gravitational  attraction 
and  centrifugal  force:  see  sec.  2.3.1. 


1  We  prefer  the  term  "kinematical'*  to  "dynamical"  for  two  reasons.  First, 
"dynamical  geodesy"  is  often  used  equivalently  with  "physical  geodesy, " 
and  second,  the  general  theory  of  relativity  is  kinematical  rather  than  dynamical 
in  character. 
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1.  Geodetic  Use  of  Gravitational  Gradients 


1.1.  Second  Derivatives 

Let  us  assume  that  all  second  derivatives  of  the  gravitational  potential  V 
in  an  earth-fixed  coordinate  system  xyz  have  been  measured  along  a  line  which, 
for  instance,  represents  the  flight  path.  These  second  derivatives  form  a  matrix 


/ 


y* 


\  v- 

where,  as  usual, 


V. 

V. 

*r 

X  z 

V 

V. 

y  y 

yz 

v.. 

V  , 

*  y 

z  z 

\ 


j£v 


ax2  ’ 

This  matrix  is  symmetric  because 

vy,  =  v*y  •  etc- 


V  .  ^ 
*y 


9xsy 


(1) 


etc. 


(2) 


If  all  second  derivatives  have  been  measured,  so  that  the  whole  matrix  (1)  is  known, 
then  the  first  derivatives,  the  components  of  the  vector  of  gravitational  force, 


F  =  (V,,  Vy,  Vz )  =  grad  V ,  (3) 

are  obtained  by  integration  along  the  flight  path.  For  instance, 

V,  -  <V„),  +  J  (V.,dx+ V.  dy+ vsidz)  . 

(VT)o  refers  to  an  initial  point  P0  on  the  flight  path,  and  Vx  refers  to  a  current 
point  P.  It  is  understood  that  the  flight  path  is  known  as  a  function  of  time  t: 
x  =  x(t)  , 

y  =  y(t)  ,  (4) 

z  =  z(t)  . 
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Then 


dx  =  x  dt ,  etc. 

(the  dot  denoting  differentiation  with  respect  to  tim^,  so  that  we  obtain 
v,  -  (V.).  +  |  (V„x+V„y+V„i)dt, 

v,  '  (V,)„  *  f  (Vy,x+ Vyyy+Vylz)  dt,  (5) 

P0 
P 

\  -  (V.).  +  J  (V„x+V>yy+V,yi)dt. 

Another  integration  gives  the  potential  V  itself: 

p 

v  =  v0  +  J*  (V,x+ Vyy+Vrz)  dt  .  (6) 

% 

The  initial  data  (Vx)0 ,  (Vy)0 ,  (Vt)a ,  V0  must  be  assumed  to  be  known. 

It  is  characteristic  for  this  method  that  the  data  (1)  are  needed  only  along 
a  line  (the  flight  path).  External  geodetic  information,  as  mentioned  in  the  Intro¬ 
duction,  is  required  in  this  method  to  obtain  the  flight  path  (4),  which  is  necessary 
to  perform  the  integrations  (5)  and  (6). 

1.2.  Second  Horizontal  Derivatives 
We  consider  the  quantities 

Wyy-  W„5WAaixlWJty  ,  (7) 

which  can  be  measured  with  a  torsion  balance  (Baeschlin,  1948;  Mueller,  1963). 
Here  W  denotes  the  gravity  potential,  the  sum  of  the  gravitational  potential  V 
and  the  potential  of  centrifugal  force.  The  derivatives  are  taken  with  respect  to 
a  local  coordinate  system  having  a  vertical  z-axis,  so  that  the  quantities  entering 
in  (7)  are  horizontal  derivatives. 
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As  usual,  we  define  the  anomalous  potential  T  as  the  difference  of  W  and 
the  normal  gravity  potential  U: 

T  =  W  -  U.  (8) 

Since  the  normal  values  of  (7),  and  UT  y ,  are  known  (Mueller,  1963,  p.  147), 
we  obtain  from  (7)  also 

Tyy  -T„  =  TA  and  Txy  .  (9) 

It  is  well  known  that 

C  *  q  T  (10) 


(G  representing  a  mean  value  of  gravity)  gives  the  height  anomaly  (the  separation 
between  corresponding  geopotential  and  spheropotential  surfaces),  and  that 

«  *  T„  ,  1J  =  -  q  T,  (11) 


give  the  components  of  the  deflection  of  the  vertical  (the  x-axis  pointing  northward). 
By  (9)  and  (11),  the  quantities 


*7*  "  Q  y  ’ 


(12) 


are  obtained  from  torsion  balance  measurements.  The  problem  is  to  find  the 
deflection  components  £  and  r\  and  the  height  anomaly  £  by  a  suitable  integration 
of  (12). 

This  integration  is  usually  performed  by  a  numerical  method,  the  mathe¬ 
matical  significance  of  which  is  not  immediately  evident.  Mathematically,  the 
integration  of  (x2)  poses  an  interesting  and  not  too  difficult  problem  in  partial 
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differential  equations,  which  we  shall  consider  now.  Apart  from  the  theoretical 
interest,  the  mathematical  insight  gained  in  this  way  may  also  be  of  indirect 
practical  use. 

We  shall  limit  our  considerations  to  a  restricted  part  of  the  earth's 
surface,  which  may  for  the  present  purpose  be  assumed  to  be  plane.  The 
quantities  (12)  are  assumed  to  be  given  at  each  point  of  this  plane  area,  as 
functions  of  x  and  y: 

£.r  -  ih  =  <P(x.y)  , 

0 

it  ~vy  *  y)  » 

where  we  have  put 

<p(x.y)  =  'g  T*r  *  *<x>y>  =  Q  TA  ‘  (] 

The  problem  is  to  solve  the  system  of  partial  differential  equations  (13)  for  f; 
and  f],  that  is,  to  find  functions 

€  =  ?(*.y) .  v  =  n(x»y) 

satisfying  (13). 

Conceptually  it  would  even  be  simpler  to  consider  the  system  of  partial 
differential  equations 

T*r  =  “  G<p(x,y)  , 

0 

Tk*  "Ty,  =  -G0(x*y)  • 

equivalent  to  (13)  but  containing  only  one  unknown  function  T,  but  the  solution  of 


(13)  is  easier. 
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Tbe  general  solution  of  the  first  equation  of  (13)  must  have  the  form 
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i  -  J  ^(x,y»)  dy»  +  F(x) 

7  n 


(16) 


where  F(x)  is  an  arbitrary  function  of  x.  This  is  readily  verified  by  substitution. 
For  reasons  of  convenience  we  put 

F(x)  *  A  +  B(x-x0)  +  U(x)  , 

where  A  and  B  are  arbitrary  constants  and  U(x)  satisfies  the  conditions 

U(*o)  =  lyxj  =  0  (17a) 

but  is  arbitrary  otherwise. 

Thus  we  obtain  from  (16),  adjoining  an  analogous  equation  for  rj  , 


i  =  f  <rtx>y')dy'  +  A  +  B(x-x0)  i-  U(x)  , 

yo 

X 

r)  *  f  <pfx',y)dx'  +  C  +  D(y-ye)  +  V(y)  , 


(18) 


where  V(y)  satisfies 


V(y0)  *  Vy(y0)  =  0  (17b) 

but  is  arbitrary  otherwise. 

Now  the  functions  U(x)  and  V(y)  are  determined  by  means  of  the  second 
equation  of  (13).  The  insertion  of  (18)  gives 

I*  -n,  *  J*  <A,(x,yW  “  |  <pr(x',y)dx'  + 
ye 

+  B-D  +  U,(x)  -  Vy(y)  **  y)  .  (19) 

Setting  y  =  y0  in  this  equation  and  considering  (17b)  we  obtain 

X 

u* (x)  *  J  <py (x'' y0 ) dx' +  v(x. y0)  -  b  +  d  . 

x 
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Since  Ux  =  dU/dx,  this  equation  may  be  integrated  to  give 

U(x)  =  J*  <x-x')<py(x\y0)dx'  +  J  ^(x',y0)dx'  +  (D-B)  (x-x,,)  . 

X  X 

o  o 

Here  we  have  used  the  fact  that 

~  J*  (x-x')  f(x')dx'  =  J  f(x')dx'  (20) 

«ZX  m  | 

and  considered  (17a).  In  exactly  the  same  way  (by  setting  x  ■  x„  in  (19)  and 

integrating)  we  find  V(y),  so  that  (18)  becomes 
r  * 

4  *  J  <P(x> y')dy'  +  J*  (x-x')  <p,  (x',y0)dx’  + 

»o  *0 

X 

+  J*  ^(x'or^dx'  +  A  +  IKx-xJ  , 

*o 

n  -  J  <p(x',y)dx'  +  J  (y-y'J^K.yW  - 
*0  »0 

y 

-  J  $(Xo  ,y')dy'  +  C  +  B(y-ye)  . 


To  further  investigate  this  solution,  we  form  the  derivatives 
y 

J 

y, 


ir  =  J  <P,(x,y')dy' +  J  <py(x\y0)dx'  +  tf>(x,y0)  +  D  , 

y  * 

1 o  o 

*  y 

V,  =  J  <py(x',y)dx'+  J  <px(x0,y')dy' -^.y)  +  B  . 


Now 


J  <P,(x,y')dy'  -  J  <px(\, y')dy'  =  J  J  <^x  (x'.y')dx'dy' 

Jo  *0 


*  y 


and 


(21) 


(22) 


B  -  D  =  , 


(23) 
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as  may  be  seen  by  setting  x  =  Xo ,  y  =  ya  in  (19).  Thus  we  find  from  (22) 

4, -v,  S  te*(x,.y*) -<pyr(x,.y,)]dx'dy  + 

*o  To 

+  «/»(x,  y0 )  +  0(xo ,  y)  -  Wx* ,  y0 )  . 

With  the  identity 

Wx,y0)  +0(xo.y)  -^.Yc)  ■  W*.y)  -  J*  f  ^(x'.y'jdx’dy* 

*e  T0 

this  becomes 


4*  -  Vr  ■  W*.  y)  +  J*  /  *  (x’ ,  y ' ’)  -  <P,  ,  (X',  y 1 ' )  -  4  y  (x* ,  y ' '>1  dx'dy '  . 

xo  To 

This  agrees  with  the  second  equation  of  (13)  only  if 


<Pxx  -<Pjy  5  0  «  (M> 

The  equation  (24)  thus  constitutes  the  integrability  condition  of  the  system 
(13).  The  data  functions  <p(x,  y)  and  ip(x,  y)  cannot  be  prescribed  arbitrarily  but 
must  satisfy  the  integrability  condition  (24).  Errorless  data  will  automatically 
satisfy  (24)  because  by  (14) 


Px*  -<Pyy  -  4y  “  r«  ('^«yj«  +  ^"xyyy  ”  ^yyxy  +  ^xxxy) 


=  o  (*^xxxy  +  ^"xyyy  ”  ^xryy  +  ^"xxxy)  «  ® 


With  empirical  data,  (24)  must  be  enforced  by  an  adjustment  if  necessary. 

To  determine  the  constants  A,  B,  C,  D  in  (21),  we  put  x  =  x„ ,  y  =  y0 
in  (21)  and  (22).  Using  (23)  we  readily  obtain  with  =  £(x„  ,y0),  etc. : 

A  -  4o.  B  =  (4x)o>  C  =  Vo  >  D  ■  (Vy)o  * 

Only  three  of  these  constants  are  independent,  because  B  and  D  are  connected 
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by  (23),  so  that 

D  =  (t?y) o  =  (4,) o  . 

Thus  (21)  assumes  the  final  form 

4  *  4o  f  HI,).  -0(Xo.yo)l  (x-xj  + 

r  *  * 

+  J*  (P(x,  y')  dy'  +  ,f  (x-x*)  <py(x*,  y0)  dx'  +  J*  0(x' ,  y0 )  dx' , 

*0  *0  *0 

(25) 

v  =  Vo  +  (l,)0(y-y0)  + 

X  y  y 

+  /  ^x'.yjdx' +  /  (y-y')^  (Xo.y’ldy1  +  J*  ^.y'Jdy*  . 

*0  y0  r0 

The  three  quantities  £0 ,  T)a ,  (4X  )0  must  be  given  at  the  initial  point  P0 , 
and  the  date  functions  <p  and  ip  must  satisfy  the  integrability  coriition  (24) . 

Then  4  and  r\  are  given  by  (25).  If  we  wish,  we  may  use  4  and  77  so  obtained 
to  get  £  or  T  by  the  usual  method  of  astronomical  leveling: 

p 

C  -  Co  “  J  (4dx  +  rjdy)  .  (26) 

po 

It  should  be  noted  that  the  integrations  needed  to  determine  4  by  (25)  are 
extended  over  the  horizontal  or  the  vertical  part  of  the  full  line  in  Fig.  1,  whereas 
the  integrations  for  r\  involve  the  broken  line  in  Fig.  1. 

This  method  presupposes  the  functions  <p(x,  y)  and  0(x,  y)  to  be  given  at 
every  point  of  a  certain  region,  thus  sharing  a  general  characteristic  of  physical 
geodesy.  It  is  theoretically  rigorous,  apart  from  the  approximation  of  the  reference 
ellipsoid  in  a  small  region  by  a  plane.  Extensions  to  a  spherical  or  ellipsoidal 
surface  of  reference,  although  simple  in  principle,  lead  to  awkward  formulas. 
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As  a  matter  of  fact,  measurements  of  <p  and  n>  are  performed  at  discrete 
points  only,  so  that  the  integrals  in  (25)  must  be  replaced  by  sums.  It  seems 
preferable,  however,  to  start  from  rigorous  expressions  such  as  (25)  and  find 
practical  approximations  for  them,  rather  than  starting  from  approximations 
right  away. 

It  is  clear  that  instead  of  prescribing  £0 ,  rj0 ,  (£x  )0  it  is  better  in  practice 
to  prescribe  £0,  r)0  and,  say  at  a  different  point  Pt.  The  expression  of  , 

Vo .  (C,)0  by  Co »  Vo  •  Ci  or  by  ,  %  by  means  of  (25)  is  straightforward 

and  may  be  left  to  the  reader. 

Finally  we  remark  that  the  horizontal  derivatives  (9)  may  also  be  measured 
by  airborne  equipment.  If  these  quantities  are  measured  in  this  way  along  a  level 
surface,  formulas  (25)  and  (26)  may  be  applied  to  get  T;  by  adding  the  normal 
gravitational  potential  to  T,  the  potential  V  is  obtained. 

1.3.  Second  Vertical  Derivatives 

Using  the  notations  of  the  preceding  section,  we  consider  the  quantity 
T 

x* t  • 

Since  the  z-axis  coincides  with  the  local  vertical,  this  is  the  second  vertical 
derivative  of  the  anomalous  potential  (8). 

Assuming  Tt  t  to  be  known  at  every  point  of  a  level  surface,  we  shall 
now  derive  an  integral  formula  expressing  T  in  terms  of  TXI .  We  consider  the 
spherical  approximation,  that  is,  we  replace  the  reference  ellipsoid  formally  by 
a  sphere;  this  approximation  is  also  used  in  Stokes'  formula. 
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Then  we  have 


T 


Z  Z 


32t  _  afr 

az2  =  ar2  "  ’ 


(27) 


since  the  normal  to  a  sphere  is  directed  along  the  radius  vector  r.  We  shall  now 
follow  a  procedure  described  in  (Heiskanen  and  Moritz,  1967,  pp.  88  and  97),  us  mg 
the  same  notations. 

We  expand  the  anomalous  potential  T  into  a  series  of  spherical  harmonics: 
T(r,  0,X)  =  £  (^T1  T„  (0,  X)  •  (28) 

n=  0  ^  r ^ 

Differentiating  twice  with  respect  to  r  we  obtain  in  space 
Trr(r,fl,X)  =  £  (n+1)  (n+2)  Tn  (0,X)  , 

0  ” 

and  at  the  level  surface  under  consideration  (r-=R) 


Trr  =  ^  £  (n+1)  (n+2)  Tn  (0,X)  . 


(29) 


We  may  also  directly  egress  Tpp  as  a  series  of  Laplace’s  surface  harmonics: 
TPP  =  £  Trp  ,(0,X)  . 

0  ’ 

Comparing  these  two  series  yields 

Trr  n  (6»X)  =  (n+1L(n+2)  Tn  (fl,X)  ,  Tn  RS 


(n+1)  (n+2) 


so  that  on  the  level  surface  (r=R) 

T  -  5 T- '">«  -  5  - 

According  to  a  well-known  formula  for  spherical  harmonics  (Heiskanen  and 
Moritz,  1967,  p.  30)  we  have 


(30) 
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rr.» 


zn-t-i 
4  n 


JjTrrPn  (COS^da 
a 


so  that  (30)  becomes 


2n+ 


0  » - /  \ - /  Q 


Bjr  interchanging  the  order  of  summation  and  integration  (which  may  be  shown  to 
be  permissible)  we  get 


2n+l 


o  (n+1)  (n+2) 


Pa  (cos$)  ]  Trrda  . 


On  putting 


s.«>  *  *  p-  «-♦>  <3i> 

we  have 

T  -  5  I!  TrA  <♦>  <b  .  (32) 

a 

This  integral  formula  is  completely  analogous  to  Stokes'  formula:  it 
expresses  T  in  terms  of  Tpr  just  as  Stokes'  formula  expresses  T  in  terms  of 
Ag.  Even  the  expression  (31)  for  the  function  S1(^i)  is  comparable  to  the  expression 
of  Stokes'  function  S(^i)  in  terms  of  spherical  harmonics  (Heiskanen  and  Moritz,  1967, 


p.  97). 

We  shall  now  develop  a  closed  expression  for  the  function  Sx(^)  by  summing 
the  series  (31).  The  function 

fM)  =  7  =  Vl-2tcos0+t2  *33* 

may  be  expanded  as 


L6 


(34) 


f(M)  ■  t  1?PB  (cos^j) 

0 

(Heiskanen  and  Moritz,  1967,  p.  33).  By  consecutive  integrations  with  respect  to 
t  we  obtain 


MM)  55  /  *M)A  =  £  p»  <C0SW  * 


4(x,0)  =  (MM>*  -L  Pn  . 


MM)  =  f  MM)dt  =  E 


P.  . 


j  -  -  (n+1)  (n+2)  (n+3)  '■  ‘ 

According  to  (Smirnow,  1966,  p.  46)  we  have 

MM)  =  J  M  f  (t)dt  , 

0 

MM>  *  £  |  (x-t)sf(t)dt  . 

From  (35)  and  (36)  we  obtain  with  x  =  1  and  l/l  defined  by  (33): 

wtsa-lM!*’ 

t  /l4  •  P,(cos^)  1  n  dt  *  1  .  1  p  &  ^ 

h3^)  =  J  (n+l)  (n+2)  (n+3)  =  *  /  X  "  /  A  *  /  X*1 

The  integrations,  to  be  performed  by  standard  methods,  give 

1 \(i>)  =  xn(l+l/sin|)  , 

hgty)  =  2  sin3 1  xn(l+l/sin|)  +  1  -  2sin|  , 

hgto)  =  (-sin3 1  +  3siii*|)  xn(in/sin j|)  + 

+  i +  f  sinS2  “  3sin3l  • 
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After  these  mathematical  preliminaries  we  can  at  once  determine  the  sum  (31). 

We  have 

■>**  = 

"  SCn+l  "  (n+1)  (n+2)  ") P"  <cos® 

=  2\Q>)  -  311,(0)  . 

so  that  from  (38)  we  obtain  our  final  result 

8.(0)  =  (2-6sin2£)  *n  (1+l/sinf)  -  3  +  68inJ  .  (o9) 

2  2  2 

This  function  is  to  be  used  in  the  integral  formula  (32) . 
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with 


os  Oq+I 

=  “  J  (n+l)  (n+2)  (n+3)  P"  (C0S^ 
corresponding  to  (31) . 

The  series  (43)  may  be  readily  summed.  We  have 

Sgty)  =  L  (“  (n+i)  (nr 2)  +  (n+l)  (n+2)  (n+3))  P»<cos^ 

-  -  2hj(^i)  +  5hg(^)  , 

using  the  notation  of  (37).  From  (38)  we  obtain  at  once  our  final  result 
S2(0)  -  (~9sin3|  +  lSsid4!)  ^(l+l/sinl)  - 

-  f  +  4sin^  +^sin2|  -  15sinp|  . 

This  function  is  to  be  used  in  the  integral  formula  (42). 


(43) 


(44) 


1.5.  Transformation  of  Gradients 

In  conformity  with  geodetic  usage,  we  shall  use  the  term  "(generalized) 
gradients"  for  the  second  and  third  derivatives  of  the  potential.  The  problem  of 
transforming  these  gradients  from  one  rectangular  coordinate  system  to  another 
occurs  frequently;  for  instance,  it  may  be  required  to  transform  them  from  a 
system  fixed  to  the  measuring  equipment  to  a  system  fixed  to  the  earth. 

It  is  convenient  to  use  the  symbols  x1 ,  Xg ,  Xg  for  the  coordinates 
x,  y,  z.  Let  the  original  system  be  (Xj\  Xg',  Xg')  ,  and  the  new  system 

(Xj ,  Xj.  Xq).  The  transformation  equation  may  be  concisely  written  as 

3 

x,’  =  £  aj  >  Xj  ,  (45) 

Jrl 
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or  as  the  matrix  equation 


h') 

l*ii 

**18 

\ 

l<\ 

1  \ 

»<? 

= 

Si 

*33 

*33  ' 

j  x= 

\  ®3l 

*33 

^33  i 

1*3  / 

Here  the  element  a{  is  the  cosine  of  the  angle  between  the  old  xt'-axis  and  the 
new  Xj-axis. 

The  matrix  (at  3 )  is  orthogonal,  that  is,  its  elements  satisfy  the  six 
conditions 


(45f) 


Lalkau 


3  *  i  . 
j  i  i  • 


(46) 


Since  the  nine  elements  a,  3  are  connected  by  six  conditions,  they  must  be  functions 

of  three  independent  parameters,  corresponding  to  the  three  degrees  of  freedom  of 

© 

rotation.  For  these  parameters  one  may  take  the  three  Eulerian  angles  6,  <p, if) 
in  terms  of  which  we  have  (e.g. ,  Goldstein,  1950,  sec.  4-4). 
a^  =  cos^cosjfc  -  sinpsin^cosd, 
a12  =  sin<pcos$  +  cos  <p  sin  ifi  cos  0, 
a^  *  sin^sinb, 

aSl  =  -  cos<psin ip  -  sin<pcos(()  cos0,  (47) 

Etgg  =  -  sincpsin^  +  cospcos(i)cos0, 
as3  =  cos^sin©, 

a3l  =  sin<psin0,  agg  =  -  cos<psin0,  ag3  =  cos 0  . 

The  gradients  in  the  original  system  will  be  denoted  by  a  prime,  so  that 
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(48a) 


oxjaxj 

similarly  in  the  new  system 


ax,'  sxj  axj 


1J  "  gXjBX,  ’  nk  r)Xj  BXj  ?xk  ’  1 

By  the  usual  rules  of  partial  differentiation  we  have 

|t,  i  jt  »£  _  _s£ 

9*t  3*J  8*1  t  ”  8Xp' 

using  (45).  Further  differentiation  yields,  on  introducing  the  abbreviations  (48a,  b), 


T,  j  «  2  aqj  Tp<i  » 

p,q=l 


T1JV  “  1*  Sia,jarVTp,r  *  <50> 

P,q,r=l 

As  an  example,  in  sec.  1.3  we  used  Tzz  ,  and  in  sec.  1.4  we  used  Tzzz. 
Since  in  our  new  notation  z  =  Xg ,  we  have  by  (49)  and  (50) 


-  T33 


ap3aq3,^pq  > 

p,  q=i 


Tr* i  -  t333 


^  ap  3  aq3  ar  3  "^p  qr  » 

P,  1,  r=1 


where  a^g  is  the  cosine  of  the  angle  between  the  old  x^  -axis  and  the  new  Xj-axis, 
that  is,  the  cosine  of  the  angle  between  the  local  vertical  and  the  old  xj  -axis. 
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2.  Kinematical  Measurement  of  Gravitational  Gradients 
(Theoretical  Background) 

2.1.  Separation  of  Gravitational  and  Inertial  Effects 

A  gravimeter  mounted  aboard  an  airplane  will  measure  the  resultant  of 
gravitational  force  and  inertial  acceleration  due  to  irregular  motion  of  the  air¬ 
craft;  the  gravimeter  will  not  be  able  to  separate  gravitational  and  inertial  com¬ 
ponents.  The  reason  is  the  equality  of  "gravitational  mass"  and  "inertial  mass, " 
which  has  been  established  experimentally  to  an  extremely  high  accuracy  by  Etttvbs 
and  others  and  explained  by  Einstein  in  his  general  theory  of  relativity.  This  is 
the  Principle  of  Equivalence,  according  to  which  gravitation  and  inertia  are 
essentially  identical. 

This  principle  seems  to  exclude  any  hope  of  separating  gravitational  and 
inertial  effects.  In  fact,  if  the  force  vector  at  one  particular  point  only  is  con¬ 
sidered,  this  force  cannot  be  separated  into  its  gravitational  and  its  inertial  parts 
because  of  the  identity  of  gravitational  and  inertial  mass.  Such  a  separation  is, 
however,  possible  with  higher-order  gradients  because  of  the  structural  differences 
between  gravitational  and  inertial  fields. 

The  practically  simplest  approach  to  this  problem  is  in  terms  of  classical 
mechanics  (sec.  2.2  ),  but  this  approach  is  not  quite  satisfactory  because  gravita¬ 
tional  and  inertial  forces  are  treated  as  different  concepts  from  the  outset,  /hich 
seems  to  be  contrary  to  the  Principle  of  Equivalence. 

A  deeper  understanding  requires  a  relativistic  analysis  (sec.  2.3  ).  The 
general  theory  of  relativity  is  considered  the  best  theory  of  gravitation  available 
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at  present.  It  gives  an  explanation  of  the  equality  of  gravitational  and  inertial 
mass  by  showing  that  these  two  kinds  of  mass  are  identical,  whereas  according 
to  classical  mechanics  this  equality  is  an  accidental  coincidence.  Although  the 
relativistic  treatment  will  confirm  the  results  of  the  classical  mechanics,  it  is 
necessary  in  order  to  dispel  the  doubts  raised  by  the  Principle  of  Equivalence. 

In  addition,  the  relativistic  approach  is  quite  appropriate  to  kinematical  gravimetry 
because  time,  the  kinematical  element,  is  woven  naturally  into  1  unified  space- 
time  structure. 

2. 2.  The  Approach  Through  Classical  Mechanics 

It  will  be  convenient  to  use  a  similar  notation  as  in  sec.  1.5,  denoting  the 
coordinates  x,  y,  z  by  Xj,  xz,  x3,  or  briefly  by  x, ,  where  i  ■  1,  2,  3. 

Consider  an  inertial  coordinate  system  Xt ,  that  is,  a  coordinate  system  in 
which  the  Newtonian  equations  of  motion 

Xt  =  Ft  (51) 

hold.  A  dot  denotes  differentiation  with  respect  to  time;  therefore,  X,  will  be  the 
components  of  the  vector  of  acceleration;  Ft  are  the  components  of  the  force  acting 
on  a  unit  mass.  Hence  (51)  is  the  well-known  law:  the  force  equals  the  product  of 
mass  (taken  as  unity)  and  acceleration. 

Consider  now  another  rectangular  coordinate  system  x,  rigidly  connected 
with  the  measuring  apparatus.  The  systems  x,  =  (x,y,  z)  and  X,  =  (X,  Y,  Z) 
will  be  connected  by  an  orthogonal  coordinate  transformation  of  the  well-known  form 

3 

Xj  =  E  j  Xj  +  Bj  ,  (52) 
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which  differs  from  (45)  only  by  B,  which  denotes  a  displacement  of  the  origin  of  the 
system  x,  corresponding  to  the  motion  of  the  aircraft.  Both  the  rotation  matrix 
a,  and  the  displacement  vector  Bt  are  functions  of  the  time  t. 

We  shall  now  introduce  the  familiar  summation  convention  (e.  g. ,  James 
and  James,  1959,  p.  377),  according  to  which  the  mere  repetition  of  an  index  is 
sufficient  to  denote  summation  with  respect  to  this  index  over  its  range.  Hence 
(52)  may  be  abbreviated  as 

Xj  8  j  Xj  +  B(  ,  (52*) 

the  summation  being  indicated  by  the  repetition  of  the  index  j. 

Repeated  differentiation  of  (52')  with  respect  to  the  time  t  yields 

•  • 

=  ^  .1  XJ  +  ai  J  X3  +  * 

(53) 

*  El  3  ^3  +  ^a»  3  *3  +  1  X3  +  * 

Combining  (51)  and  (53)  we  obtain 

ai3x3  *  Fi  ”  2ai  j *3  “a«jx3  -B*  • 

Multiplying  by  a,  k  (the  summation  with  respect  to  i  is  then  automatically  implied) 
we  obtain 

ai  k  ai  3  X3  =  ai  k  FJ  k  ai  1  X3  “  af  k  "at  j  XJ  ”  ai  k  ^ 

In  agreement  with  (46)  we  have 

aika!3  =  6V3  (55> 

(the  summation  may  as  well  be  over  the  first  index  i  as  over  the  second  index  as 
in  (46)).  Hence  the  left-hand  side  of  (54)  becomes 

ai  k  ai  3  X3  =  ®k  3  X3  "  \  ’  (56a> 

since  6k }  is  zero  for  j  4  k  and  unity  for  j  =  k. 


Multiplying  (52’)  by  alk  we  find 
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The  matrix  Wj  k  characterizes  the  instantaneous  rotation  since  it  may  be  shown  that 
the  vector 

u)  =  (coj,  ufe,  w3)  ; 

COj  “  Wgg  ,  OJg  —  W3,  ,  (O3  ~  Wjg 

has  the  direction  of  the  instantaneous  axis  of  rotation  (with  respect  to  the  moving 
system  x, ) ;  the  magnitude  of  oj  , 

|  W  |  =  ^vfa  +  wfi  +  wfa  =  to  , 
is  the  instantaneous  angular  velocity  between  the  systems  x,  and  Xj  . 

By  differentiating  (60)  we  find 

y  ~  v  ai  s  +  at  v  j  ; 
on  the  other  hand, 

WJjtWk  t  ~  ~  \nKi\v  =  a«ja!k  • 

so  that 

aivaU  =  WJV  -  wUwvX  •  <62> 

With  (60)  and  (62),  equation  (58)  goes  over  into 

K  =  K~  -(wJk  -Bk 

or,  on  replacing  the  index  k  by  i  and  using  (61), 

x,  =  f,  +  2w, }  Xj  +(wi}  +wfkwJk)x3  -B,  .  (63) 

This  is  the  final  form  of  the  equations  of  motion  in  the  moving  system  x,  -1 
It  is  not  of  the  Newtonian  form  (51).  To  put  it  into  a  Newtonian  form, 

x,  =  f,*  ,  (64) 

For  a  more  sophisticated  derivation  see  (Morgenstem  and  Szab6,  1961,  pp.  7-9). 
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we  must  set 


f  *  =  fi  +  2Wi  3^3  +  <W1  kwjk  +  w,  3)Xj  -  bj  ,  (65) 

thus  adding  to  the  '•true”  force  f,  certain  fictitious  "inertial  forces, "  of  which 

2W1  3^  3 

is  the  Coriolis  force,  and 

WlkW^Xj 

is  the  centrifugal  force. 

If  we  consider  f,  to  be  the  gravitational  force,  it  is  the  gradient  vector  of 
the  gravitational  potential  V, 


f  =-^ 
1  ax, 


(66) 


but  the  output  of  our  measuring  system,  which  operates  according  to  the  Newtonian 
form  (64),  will  be  f,*.  By  (65)  and  (66),  with  x,  =  0  because  the  measuring  system 
is  at  rest  with  respect  to  the  frame  xyz,  we  have 


3V 

dx, 


=  f*  -(wlkwn  +wn)Xj  +  B, 


(67) 


This  equation  forms  the  starting  point  for  our  conclusions.  Successive 
partial  differentiations  with  respect  to  x5  yield 


d5V 

9x,  ax3 

a3v 


af* 

9x, 


"  <wlkwJk  +  w13)  , 
a2f* 


(68) 

(69) 


9x1^xJ^xk  aXj^Xk 

We  shall  now  distinguish  two  cases : 

1.  No  inertial  stabilization.  The  measuring  instrument  is  rigidly  connected  to  the 
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aircraft  and  follows  its  linear  acceleration  and  its  rotational  movements.  Then 


both  the  gravitational  vector,  by  (67),  and  the  second  derivatives,  by  (68),  will 
falsified  by  inertial  effects,  but  the  third  derivatives  of  the  potential  will  no 
longer  be  affected  by  inertial  disturbances,  as  (69)  shows. 

2.  Inertial  stabilization.  Here  the  instrument  axes  xyz  will  be  kept  constantly 
parallel  to  inertial  axes  XYZ  (to  assume  the  simplest  case)  by  inertial 
stabilization.  Then 
w, ,  SO. 

and  (67)  and  (68)  reduce  to 

av 


ax, 


-  f  *  +  S,  , 


a3v  , 


(70) 


(71) 


ax,  axj  axj 

Here  the  gravitational  vector,  by  (70),  is  still  affected  by  linear  accelerations, 
but  inertial  disturbances  are  absent  already  in  the  second  derivatives  (71). 

Hence,  by  measuring  higher  derivatives  of  the  potential,  it  is  indeed  possible 
to  separate  gravitational  effects  from  inertial  disturbances.1 


1  Preliminary  considerations  along  this  line  have  been  made  by  Balabushevich  (1954); 
the  separability  in  the  higher  derivatives  was  also  recognized  by  Veselov  (1964) 
through  a  different  reasoning. 
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2.3.  The  Approach  Through  the  General  Theory  of  Relativity 
2.3.1.  Principles;  the  Kinematicai  Potentials 

In  the  theory  of  relativity,  space  is  combined  with  time  into  a  space-time 
continuum.1  The  space-time  of  the  special  theory  of  relativity  is  a  four  dimensional 
manifold  with  coordinates 

x1=x,  Xs  -  y  ,  X3  =  z  ,  x4=t,  (72) 

whose  line  element  ds  may  be  expressed  as 

ds2  *  (^dt?  -  (dx3  +  dy3  +dz2)  , 

where  c  is  the  velocity  of  light  in  vacuum.  For  the  purpose  of  the  present  paper 
it  is  more  convenient  to  divide  by  the  constant  cs,  thus  using  for  the  line  element 
the  expression 

ds2  =  dt2  -  c-a(dx2  +  dy2  +  dza)  .  (73) 

This  may  also  be  written 

ds2  *  dt2  -  c”adx,dx1  ,  (74) 

where  Latin  indices  denote  the  spatial  indices  1,  2,  3  only;  by  (72)  this  expression 
is  identical  with  (73),  the  summation  convention  being  kept  in  mind. 


1  There  are  many  excellent  introductions  to  the  theory  of  relativity.  We  mention 
(Bergmann,  1942),  (Adler,  Bazin,  and  Schiffer,  1965)  and,  more  advanced  but 
highly  original,  (Synge,  1964  and  1960);  some  more  may  be  found  in  the  list  of 
references. 
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Now  we  shall  introduce  a  further  formal  change,  writing  superscripts  instead 
of  subscripts: 

ds3  =  dt3  -  c^dx*  dx1  .  (74') 

Here  dx1  is  only  another  notation  for  dx,  ;  this  notation  is  in  agreement  with 
tensor  calculus,  where  "co variant"  indices  are  distinguished  from  "contravariant" 
indices,  the  former  being  written  as  subscripts,  the  latter  as  superscripts. 

Because  the  dx1  form  a  contravariant  vector,  the  notation  (74')  is  in  agreement 
with  tensor  notation  and  therefore  sometimes  preferable. 

The  expression  (74')  has  the  general  form 

ds3  =  g^dxW*  ,  (75) 

if  we  include  dt  as  dx?  in  agreement  with  (72)  (note  that  the  superscript  4  is  not 

an  exponent!)  and  let  Greek  indices  run  from  1  to  4,  whereas,  as  we  have  mentioned, 

Latin  indices  run  from  1  to  3  only.  Because  of  the  summation  convention,  (75)  is, 

of  course,  equivalent  to 

*  4  at 3 

dsa  =  t  L  dx  . 

c*=l  p=l  ttfr 

fhe  matrix  gap  ,  the  covariant  fundamental  tensor,  has  the  form 


0  0  0  1 


as  we  recognize  by  comparing  (75)  and  (74').  That  is,  we  have 
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0 


(77) 


St4  "  1  . 

$4  "  g41  = 

g13  »  -c-s5 

The  fundamental  tensor  g ap  has  this  simple  diagonal  form  with  constant 
coefficients  only  in  an  inertial  coordinate  system  without  gravitational  forces. 

Such  inertial  coordinate  systems  in  space-time  are  the  four-dimensional  analogue 
of  rectangular  coordinates  in  space. 

If  we  use  curvilinear  coordinates  in  space-time,  then  the  g^  will  no 
longer  be  equal  to  the  constant  values  (77).  This  will  happen  with  non- inertial 
coordinate  systems  such  as  accelerated  or  rotating  frames:  although  the  spatial 
3-system  is  rectangular,  it  forms  together  with  time  a  curvilinear  4-system. 

As  an  example,  consider  the  rectangular  non- inertial  system  x,  =  (x,  y,  z)  of 
sec.  2. 2.  hi  the  inertial  system  X.  =  (X,  Y.  Z)  ,  the  fundamental  tensor  will 
have  the  simple  form  (76),  so  that 

ds2  =  dt2  -  c-2dXj  dX,  .  (78) 

According  to  the  transformation  formula  (52')  we  have 
dX,  =  a,  j  dXj  +  (a, }  x}  *  Bj )  dt  . 

Substituting  this  into  (78)  we  obtain 

ds2  =  ga/3dxadx^  (79) 

where  now 

g,4  =  l  -  cr2(wjjewkiexJxk  -  2wJkbJxk +b}b3)  , 

Sj 4  =  g*,  =  crs(wnx3  -  bj)  ,  (80) 

cr  —  —  /*-  2  6 

c  u* 
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where  k  is  the  gravitational  constant,  dm  is  the  element  of  mass,  and  i  is  the 
distance  between  dm  and  the  point  to  which  V  refers. 

It  may  be  mentioned  that  the  h4 }  are  also  related  to  V  by 

s 

ht i  =  •2V5tJ  . 
s 

This  fact  will  not  be  used,  however,  because  in  this  paper  we  shall  consistently 
limit  ourselves  to  the  linear  approximation  (linear  in  c“2)  ,  which  at  present  is 
sufficient  for  almost  all  practical  purposes: 
g44  =  1  -  2c“ 2  V  , 

*i  =  0  ,  (85) 

Si  j  =  C"261  j  • 

The  fact  that  only  the  Newtonian  potential  V  occurs  in  the  linear  approxi¬ 
mation  suggests  that  the  results  obtained  by  using  (85)  will  be  essentially  equivalent 
to  classical  mechanics;  and  every  textbook  on  general  relativity  shows  that  this  is 
indeed  true.  Physically,  the  linear  approximation  corresponds  to  a  slow  motion 
(velocity  «c)  in  a  weak  gravitational  field  (such  as  that  of  the  earth,  the  sun,  etc.). 

Let  us  now  perform  the  transition  to  an  accelerated  and  rotating  frame.  In 
exactly  the  same  way  as  we  transformed  (77)  into  (80),  we  transform  (85)  into 
g*4  =  1 -c-a(2V  +  wj£wk£xjXk  -2wlkbjXv -t-bjbj)  , 
g14  =  g»i  =  c"2(w13Xj  -b,)  ,  (86) 

*  -c"26n  * 

To  get  a  better  understanding  of  these  formulas,  consider  a  system  xyz 
connected  with  the  rotating  earth,  the  z-axis  coinciding  with  the  axis  of  rotation. 

33 


I 


Then 


b,  sO, 


<wn)  =  l 


0  to 

00  o 
0  0 


so  that  (86)  reduces  to 


»/ 


8*4  =  1  -  c"a[2V  +  uFfx2  +  y2)]  , 
gi4  ‘  <rawy  , 

Bb4  *  “  c"a  1DX  , 


(87a) 


(87b) 


(88) 


834  “  ®  • 
g,.  -  ~  c“a6n  . 

Introducing  the  potential  of  gravity  (gravitation  plus  centrifugal  force), 

W  -  V  +  ioo2(xa  +f)  ,  (89) 


we  see  at  once  that 


g*4  -  1  -  2c“a  W  ;  (90) 

that  is,  gj4  is  (apart  from  constants)  essentially  identical  with  the  potential  of 
gravity. 

Similarly,  if  there  is  no  rotation,  (85)  shows  that  g^4  is  essentially  identical 
with  the  potential  of  gravitation,  V.  The  fact  that  in  (88)  the  potential  of  centrifugal 
force  enters  into  g^4  in  a  natural  way,  is  another  expression  of  the  Principle  of 
Equivalence,  since  gravitational  and  inertial  forces  occur  in  g^4  side  by  side  in 
an  equivalent  manner. 
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Thus  in  (88),  contains  the  statical  components  of  the  gravity  field: 
gravitation  and  centrifugal  force.  It  may  be  shown  that  the  velocity-dependent 
component  of  the  gravity  field,  the  Coriolis  force,  is  expressed  by  the  quantities 
g,4 .  We  therefore  infer  that,  also  in  the  general  case  of  (86),  the  complete  set 
of  gap  will  express  all  features  of  the  combined  gravitational-inertial  field;  the 
gap  may  therefore  be  called  the  kinematical  gravity  potentials. 

Thus  the  extension  of  physical  geodesy  to  moving  systems  requires  the 
introduction  of  the  complete  set  of  the  kinematical  potentials  gap,  the  component 
g^4  equivalent  to  the  statical  potential  W  is  no  longer  sufficient  to  characterize 
the  gravity  field. 

That  the  set  of  is  necessary  and  sufficient  for  this  purpose  is  expressed 
by  the  fact  that  the  path  of  a  particle  under  the  combined  influence  of  gravitational 
and  inertial  forces  is  a  geodesic 


J ’  * 


minimum  . 


(91) 


The  path  is  therefore  determined  by  the  gap  .  We  shall  return  to  this  point  later, 
but  we  mention  already  here  that  the  evaluation  of  (91)  by  means  of  the  calculus  of 
variations,  using  the  ga ^  as  expressed  by  (86),  leads  just  to  the  equations  of  motion 
(63)  found  previously  (the  restriction  to  the  linear  approximation  being  always  under¬ 
stood). 

Velocity,  acceleration,  force.  Since  space-time  is  four-dimensional,  it 
is  natural  to  define  four-dimensional  equivalents  of  terms  such  as  velocity, 
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acceleration,  or  force.  The  4-velocity  is  simply  defined  as 

a  _  6xa  S  dx  dy  dz  dt  x 
U  ds  =  ds  *  ds  *  ds  ’  ds  ) 

From  expressions  such  as  (78),  (80),  or  (82)  we  see  that  in  every  case 

ds2  =  dt?  +  0(crs)  , 

where  0(cra)  denotes  small  terms  of  the  order  of  c~3.  Hence  also 


so  that 


ds  =  dt  +  0(<r3)  , 


fs-  1  + "<■>■*)  • 


Neglecting  a  relative  error  of  0(<ra)  in  agreement  with  our  linear  approximation, 


we  simply  have 


(x,  y,  z,  1)  +  0(<ra)  , 


where,  as  usual,  x  =  dx/dt,  etc.  Hence,  to  the  approximation  considered,  the 

ot 

spatial  components  of  the  4-velocity  u  are  the  components  of  the  ordinary 


3 -velocity 


v1  *  x,  =  (x,  y,  z)  ; 


the  temporal  component  is 


U4  =  1  . 


The  four-dimensional  generalization  of  acceleration  is  less  trivial.  It 

would  be  obvious  to  define  4-acceleration  by 

du"  „  d3xtt 
ds  ds3  ’ 


Cc 

in  analogy  to  (92),  but  it  is  shown  in  general  tensor  calculus  that  du  /ds  is  not  a 
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vector  (as  far  as  transformation  properties  are  concerned).  Instead,  it  is  necessary 

to  define  the  vector  of  4-acceleration  by 

.  a  „s  a 

Q!  fill  6  X  /OfiN 

a  =  —  =  TT  ,  (9b) 

6s  6s3 

where  6ua/6s  is  the  "co variant  derivative" 

6ua  dua  p  y 

6s  ds  fiy  ’  '  ’ 

the  Christoffel  symbols  being  given  by 


r«  -  i  +  *h±  _ 


'  ax'  ax^  axu  '  » 

where  the  matrix  (g“^),  the  contravariant  fundamental  tensor,  is  inverse  to  the 
matrix  (gff/3)  . 

The  second  term  on  the  right-hand  side  of  (97)  is  of  decisive  importance, 
because  it  expresses  the  gravitational  and  inertial  forces,  hi  fact,  we  have  (always 
in  the  linear  approximation),  for  the  inertial  system,  the  ga^  defined  by  (77): 
a1  -  x,  , 

(99al 

a4  *  0  ; 

for  the  purely  gravitational  system  (85)  : 

•:  9Y  , 


a4  -  0  ; 

for  the  non-inertial  gravitational  system  (86) : 


4  ** 

a1  =  x, 


-  —  -<wM  +wlkw}k)Xj  +  b,  , 


a4  =  0  ; 


(99c) 
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Comparing  this  with  (67)  we  see  that 

f1  =  -  v  •>  (104> 

that  is,  apart  from  the  sign,  the  spatial  components  of  the  4-force  f**  are  identical 
with  the  apparent  force  f*  of  sec.  2. 2.  The  difference  in  sign  expresses  the 
principle  of  action  and  reaction:  f  *  is  defined  as  the  acting  force,  whereas  f 
is  the  force  of  reaction  arising  from  the  measuring  system  being  forced  to  remain 
at  rest  in  the  frame  xyz. 

The  agreement  between  classical  and  relativistic  analysis,  as  expressed 
in  (104),  was  to  be  expected  for  the  linear  approximation,  because  Einstein's 
theory  goes  over  into  Newton's  theory  for  weak  fields  and  low  velocities.  To  a 
certain  extent,  this  fact  confirms  the  results  obtained  at  the  end  of  sec.  2. 2  , 
because  with  (104)  also  the  validity  of  (69)  is  verified.  As  for  the  case  of  inertial 
stabilization,  leading  to  (71),  the  possibility  of  such  a  stabilization  in  the  presence 
of  a  gravity  field  ought  to  be  investigated  by  means  of  the  general  theory  of  relativity. 
This  will  be  done  in  sec.  2. 3.3  ;  an  inherently  relativistic  method  will  then  provide 
a  final  justification  in  sec.  2.3.4. 


2.3.2.  Gravitational  Gradients  and  the  Riemannian  Tensor 

Let  us  start  with  a  simple  example  familiar  to  the  geodesist.  Consider 
rectangular  coordinates  x,  y  in  the  plane;  then  the  line  element  will  be  given  by 

ds2  =  dx8  +  dy2  .  (105) 

On  introducing  curvilinear  coordinates  u,v  in  the  plane  by 
x  =  x(u,  v)  , 


(106) 


y  *  y(u,v)  , 
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the  line  element  (105)  takes  the  form 

ds3  =  Edit8  +  2Fdudv  +  Gdv2  ,  (107) 

where  E,  F,  G  are  functions  of  u  and  v. 

The  line  element  on  an  arbitrary  curved  surface  has  the  same  form  (107). 

We  may  therefore  ask  by  which  criterion  one  is  able  to  recognize  whether  a  line 
element  of  the  form  (107)  refers  to  a  curved  surface  or  to  a  plane.  The  answer 
(due  to  Gauss)  is  as  follows:  Form  a  certain  expression  containing  E,  F,  G  and 
its  first  and  second  derivatives;  this  expression  is  called  Gaussian  curvature.  K 
the  Gaussian  curvature  is  zero,  then  (107)  refers  to  a  plane  and  can  be  transformed 
into  the  simple  form  (105);  if  tbe  Gaussian  curvature  is  nonzero,  then  (107)  refers 
to  a  curved  surface  and  cannot  be  transformed  (over  a  finite  area)  into  (105). 

The  remarkable  fact  is  that  the  principle  of  Einstein's  theory  of  gravitation 
is  the  extension  of  two-dimensional  Gaussian  surface  theory  to  four-dimensional 
space-time;  this  fact  should  be  particularly  appealing  to  the  geodesist. 

The  flat  space-time  of  special  relativity  corresponds  to  the  plane;  in  fact, 
by  substituting 

ic-1x  *  4  ,  ic-1y  =  tj  ,  ic_1z  =  £  , 
the  line  element  (73)  goes  over  into 

dsa  =  d$a  +  drf  +  d£a  +  dt?  ,  (108) 

which  is  the  precise  four-dimensional  equivalent  of  (105).  Thus  the  inertial 
systems  of  special  relativity  correspond  to  the  rectangular  coordinate  systems 
in  the  plane.  Transformations  to  accelerated  and  rotating  frames  after  the  fashion 
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of  (52*)  correspond  to  the  transition  to  curvilinear  coordinates  in  the  plane 
(equations  (106)),  and  the  resulting  expression  (79)  is  the  four-dimensional 
equivalent  of  (107);  in  fact,  in  two  dimensions, 

Sn  =  E»  Sis  =  fei  =  F»  *  G  .  (109) 

A  transition  to  "curvilinear  coordinates"  (accelerated  and  rotating  frames) 
in  space-time  without  gravitation  does  not  change  its  flat  character,  in  the  same 
way  as  the  plane,  even  when  expressed  in  curvilinear  coordinates,  still  remains 
a  plane.  Thus  the  gap  of  (80)  still  refer  to  flat  space,  although  inertial  forces 
occur. 

The  gap  of  (82),  however,  which  express  a  true  gravitational  field,  are  of 
a  completely  different  character.  They  are  comparable  to  the  E,  F,  G  of  a  curved 
surface.  Gravitation  corresponds  to  a  curved  space-time. 

hi  surface  theory  the  criterion  of  flatness  is  the  Gaussian  curvature.  Its 
four-dimensional  generalization  is  the  Riemannian  curvature  tensor,  a  set  of  20 
independent  quantities.  If  all  these  quantities  are  zero,  if 

Ra/3y6  =  0  »  (HO) 

there  is  no  gravitation,  no  matter  how  complicated  the  gff^  look;  in  this  case 
there  are  only  inertial  forces. 

Gravitation  gives  rise  to  a  curved  space-time  and  thus  to  a  nonzero 
Riemannian  tensor;  we  may  write  symbolically  (Synge,  1960,  p.  109)  : 

Ra0y6  "  graviteti01121!  field.  (Ill) 

The  analytical  expression  of  Ra^  ^  contains  the  gap  and  its  first  and 
second  derivatives;  we  have 


i 
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,  /  °3Y 

Raft'6  =  *  A* 


“2-  _ 


ax^ax6 


IfeiLA 

^  ol  y  / 

ax  ax 


*V(Vfr-W)’ 

the  Christoffel  symbols  being  defined  by  (98)  . 


(112) 


Since  contains  second  derivatives  of  g^,  we  cannot  hope  to  be 

able  to  separate  gravitational  and  inertial  effects  using  the  potential  (~g44)  or 
its  first  derivatives,  that  is,  the  force  vector,  only.  Separation,  if  possible  at 
all,  can  only  occur  with  second  or  higher  derivatives;  this  is  in  agreement  with 
the  results  of  sec.  2. 2  obtained  from  classical  mechanics.  This  is  another 
confirmation  of  the  importance  of  second  (and  possibly  higher)  derivatives  in 
kinematical  geodesy. 

The  individual  components  of  the  Riemannian  tensor  show  a  clear 
distinction  as  to  order  of  magnitude,  the  differences  being  of  the  order  of 
cf3 .  Using  (82)  with  (112)  we  find 

*■«♦<•  *  -  0(<rS»  • 

l  *  °!e-2>  •  U’2> 

C3®!  Jk4  =  0  (C  3)  ; 

the  other  components  are  either  expressible  through  these  or  are  zero  (e.g., 

Synge,  1960,  p.  16).  The  Latin  subscripts  denote  spatial  indices,  as  usual. 

We  see  that  R,  s „  f  and  R,  Jk4  are  smaller  than  R,  4 } 4  by  the  factor 
c”3  .  Hence,  neglecting  0(c"3)  as  usual,  we  have 
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l 

c8^ljV4  =  ®  » 

so  that  the  only  surviving  equation  of  (113)  is  the  first;  with  (83)  we  have 

i2y 


C2 434 


_ 0_ 


ax, 


(114) 


This  simple  equation  shows  that  the  component  R,4J4  of  the  Riemannian  curvature 
tensor  is  essentially  equivalent  to  a  second  derivative  of  the  gravitational  potential 
V.  This  will  be  extremely  important  for  the  following  sections. 

If  we  had  used  the  gap  of  (86)  instead  of  (85)  we  would  have  obtained  the 
same  result  (114)  since  the  contribution  of  the  inertial  forces  to  the  Riemannian 
tensor  is  zero. 


2.3.3.  Gyroscopic  Stabilization  and  Fermi  Propagation 

At  the  end  of  sec.  2. 2  we  have  seen  that  gravitational  effects  are  separated 
from  inertial  disturbances  already  in  the  second  derivatives  if  the  instrumental  axeB 
xyz  are  kept  constantly  parallel  to  the  inertial  axes  XYZ. 

This  result  is  consistent  with  the  relativistic  approach  of  sec.  2.3.1  ,  but 
the  method  of  keeping  the  instrumental  axes  permanently  parallel  must  be 
examined.  The  direction  of  the  instrumental  axes  is  kept  constant  by  ineitial 
stabilization  by  means  of  gyroscopes.  The  principle  is  that  the  axis  of  a  freely 
spinning  gyroscope  maintains  its  direction  even  if  its  frame  is  accelerated  or 
rotated;  furthermore,  the  axis  is  unaffected  by  gravity. 

It  is  by  no  means  obvious  that  this  is  true  also  in  the  general  theory  of 
relativity;  one  should  even  be  inclined  to  suspect  the  contrary.  Therefore,  this 
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problem  deserves  closer  examination. 


First  of  all,  the  question  arises  as  to  whether  something  like  constancy 
of  spatial  direction,  or  briefly,  a  nonrotating  frame,  can  be  defined  within  general 
relativity.  Surprisingly  enough,  this  is  possible  even  in  a  local  manner  (an 
example  to  the  contrary:  uniform  motion  of  a  mass  point  cannot  be  defined  in 
an  intrinsic  local  way).  The  relevant  concept  is  Fermi  propagation,  which  is 
considered  in  detail  and  used  extensively  in  (Synge,  I960),  hi  view  of  the  importance 
of  Fermi  propagation,  it  is  surprising  that  this  concept  is  not  even  mentioned  in 
any  other  standard  text  on  relativity  known  to  the  author.1 

Synge  interpreted  Fermi  propagation  physically  by  means  of  an  idealized 
optical  experiment  ("the  bouncing  photon").  Pirani  (1956)  pointed  out  that  the 
axis  of  an  idealized  gyroscope  ("spinning  test  particle")  also  follows  Fermi 
propagation.  To  show  this,  he  used  Papapetroifs  (1951)  equation  of  motion  of 
the  spin  of  a  spinning  test  particle.  Papapetrou’s  (1951)  developments  refer  to 
a  free  test  particle;  he  later  (Papapetrou  and  Urich,  1955)  generalized  them  to  a 
spinning  particle  moving  in  an  electromagnetic  field  and  found  the  same  spin 
equations.  For  our  present  purpose,  where  we  have  a  free  gyroscope  whose 
frame  performs  a  constrained  arbitrary  motion,  the  case  of  a  freely  spinning 
test  particle  under  the  influence  of  external  forces  would  be  relevant.  This  case 
can  be  handled  by  an  extension  of  a  method  of  (Weber,  1964).  The  result  is  that 

1  The  "successive  Lorentz  transformations"  of  (Miller,  1952,  46  and  96)  are 

essentially  equivalent  to  Fermi  transport  specialized  to  space-time. 
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Consider  now  a  system  of  three  orthogonal  vectors  Xp  ,  each  of  them 
being  represented  by  the  axis  of  a  freely  spinning  gyroscope.  In  this  way  the 
axes  of  a  rectangular  xyz  system  in  space  may  be  realized  physically.  Will 
now  the  respective  coordinate  axes  so  defined  remain  parallel  to  each  other  in 
the  ordinary  sense  if  the  frame  of  the  three  gyros  is  moved  arbitrarily  in  space, 
corresponding  to  the  motion  of  the  carrier  aircraft? 

The  answer  will  be  obtained  by  applying  the  linear  approximation  (85) 
to  Fermi  propagation.  For  the  spatial  components  X1 ,  in  which  we  are  interested, 
we  may  write  (118)  as 
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6s 


flu1 


(119) 


where  the  factor  $  is  given  by 


fi 
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ga/3*  6s  * 


(120) 


In  the  global  coordinate  system  xyzt  of  (85)  the  series  expansion  of  the  vector  X 
in  the  fashion  of  (82)  reads 
X1 


a 


X4  = 


X1  +  C'a  X1  +  c“4  X1  +  •  •  •  i 

o  i  a 

c"8  X4  +  c'4X4  +  •"  , 
l  a 


(121) 


because  a  non-zero  X4  would  be  incompatible  with  the  orthogonality  condition 
o 

(117).  We  further  have,  in  agreement  with  definition  (96), 
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a1  =  a1  +  c“8  a1  + 


(122) 
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because  a4  =  0  from  (99b)  . 
o 

We  now  insert  (85),  (121),  and  (122)  into  (120).  The  result  is 

*=  (1  -  2c-aV)c-4X4a4-c“2X,a,  +  C(cr4)  =  0(cra)  . 
11  0  0 


Thus  from  (119)  we  have 


ftX*  n/  -S\ 
—  =  «(o  »)  ; 


(123) 


that  is,  the  change  of  X1  (the  gyro  axis)  is  of  the  order  of  c“2  and  may  therefore 
be  safely  neglected. 

Thus,  general  relativity  confirms  the  fact  that  gyroscopic  stabilization  of 
instrumental  axes  is  indeed  possible;  the  effect  of  the  gravity  field  is  negligible 
in  this  case. 


2.3.4.  Application  of  Synge's  Method 

Synge  (1960,  pp.  156-158)  has  given  a  general  theory  of  a  "relativistically 
valid  geodetic  survey”  to  determine  the  Riemann  tensor,  which  expresses  the 
gravitational  field.  His  formula  is 


The  subscripts  (a) ,  (b),  (c),  (d)  assume  the  values  1,  2,  3,  4  each.  The  coordinates 
ot 

x  are  expressed  in  terms  of  four  parameters  y^  j  by  equations  of  the  form 

x  =  x  (y(i),  y(8).  y(3)»  y(*))  »  (125) 


then 


a  _  Sxa 

<•> =  sy(.) 


(126) 
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form  a  tetrad  to  which  the  components  R^tbc4)  refer. 

We  shall  now  specialize  the  parameters  in  such  a  way  that  the  base  vectors 


of  the  tetrad  become 


X“t)  =  for  i  «  1,  2,  3,  ; 

Ya  .a  a 

(*)  =  V)  =  u  * 


(127a) 


(127b) 


01  Ct  Ct  (X 

Here  X^j,  Xpj,  X^4j  are  to  form  a  system  of  four  orthogonal  unit  vectors, 

Gt 

the  first  three  X^ ,  lying  in  the  instantaneous  3-space  of  the  moving  particle,  and 
the  "time-like"  unit  vector  X(4)  is  identical  with  the  4-velocity  u  tangent  to  the 

Ot  01  OL 

world  line  of  the  moving  apparatus.  The  vectors  X  ^  j  ,  ^(2)  »  X  are  propagated 

along  the  world  line  by  Fermi  transport.  The  are  the  components  of  the 

01 

Riemannian  curvature  tensor  in  the  local  system  of  the  four  unit  vectors  X^aj  . 

By  comparing  (126)  and  (127a)  we  see  that 


■  X?, 


a  dx 


—  -  x(4)  =  X(4)  -  u  -  —  , 

so  that  along  the  world  line  under  consideration  we  have  (the  point  s  ■  0  being 
suitably  chosen) 


and  by  (93) 


y(4)  =  t  +  0(c“3) 


(129a) 


Similarly  we  may  take 


y(1)  -  x  +  0(c"3)  , 
y(2)  *  y  +  o(c-3)  , 
y(3)  =  z  +  0(c"3)  ; 


(129b) 
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according  to  (123),  the  change  of  the  spatial  Xj,j  is  only  of  the  order  of  c“2, 
so  that  to  this  accuracy  they  may  be  identified  with  the  unit  vectors  of  three- 
dimensional  inertial  coordinates;  therefore  the  identification  (129b)  is  per¬ 
missible. 

With  the  measuring  accuracy  possible  at  present,  it  is  best  to  concen¬ 
trate  on  the  components  R(!4J4)  of  the  curvature  tensor,  where  i,  j  =  1,  2,  3  ; 
by  (113)  these  components  are  predominant  by  the  order  of  c2 . 

The  following  derivations  are  rigorous;  the  approximations  (129a,  b) 
will  be  introduced  only  at  the  final  stage.  By  (124),  (127b)  and  (128)  we  have 
along  the  world  line 


"C1'4)  ■  g“PXW  '  6^)  • 


a 


Introducing  the  4-force  f  by  (102*)  we  obtain 
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P 


We  shall  now  prove  the  relation 
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By  the  definition  of  the  covariant  derivative  we  have 
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(130) 


(130') 


(131) 
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Because  the  order  of  two  consecutive  partial  differentiations  can  be  interchanged, 
and  because  Y®yb  is  symmetric  with  respect  to  y  and  6,  these  two  expressions 
are  identical,  and  (131)  has  been  proved. 

Hence 


jfrL . 


(132) 


The  second  derivative  on  the  right-hand  side  may  be  evaluated  by  twice  applying 
the  formulas  for  Fermi  propagation;  for  convenience  we  drop  the  subscript  (j). 
First  (118)  gives 

-  g  X8  ufl 

6  s  8y6A  6s 


or,  by  (102‘)  , 


• 

6s  /j  v 


(133) 


P 


Writing  (115)  for  an  arbitrary  vector  L  we  have 

B 

Substituting  for  L  the  expression  (133)  gives 

A'(fV-fV)  , 

6s  V  6s  J  6s2  sy6  6uv 

again  using  (102').  Since  u^  is  a  time-like  unit  vector  we  have  (Synge,  19C0, 
p.  2) 


g  uyufi  =  -  1 

v6 


and  on  differentiation  with  (102') 


g  .fyu6  *  0  . 
y  6 


(134a) 


(134b) 
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By  means  of  these  equations  we  finally  obtain 
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(135) 


Hence  (132)  becomes 
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uv 
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so  that  (130’)  goes  over  into 

V»>  *  *„/<•>  6^  '*«/(•) ^uVX%&-  <136> 

We  now  have  (Synge,  I960,  p.  10) 

Vw  ■  *  fo  •  (137) 

the  "covariant"  components  of  f°  with  respect  to  the  orthonormal  triad  X^  . 

By  differentiation  of  (137)  we  find 
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The  first  term  on  the  right-hand  side  is  zero  because  the  covariant  derivative 
of  the  fundamental  tensor  g  is  identically  zero;  the  second  term  can  be  made 

ur 

zero  by  a  suitable  choice  of  parameters  y^  ,  and  the  third  term  is  identical 
with  the  first  term  on  the  right-hand  side  of  (136).  Thus  (136)  becomes, (137) 


also  being  considered, 


This  simple  result  is  rigorous.  Still  using  no  approximation,  it  is 
appropriate  to  go  over  to  upper  ("contravariant")  indices  by 

f(.)  -  e(,.)J<->  , 

where  g^mj  is  the  matrix  (76);  cf.  (Synge,  1960,  p.  9).  Since 

f(4)  =  f  u"  =  0 
a 

by  (134b),  we  find 

fyj  =  -  cr2f  0)  .  (139) 

Thus  (138)  becomes 

B  “  c“af(*)f (0  .  (140) 

This  formula,  slightly  less  simple  than  (138)  and  as  rigorous,  is  more  instructive 

e 

because  f  (0  is  equivalent  to  an  ordinary  3-force  also  with  respect  to  magnitude,  and 
thus  the  orders  of  magnitude  are  immediately  obvious  in  (140). 

Now  we  shall  finally  proceed  to  our  usual  linear  approximation.  Evidently 
the  second  term  on  the  right-hand  side  of  (140)  can  be  neglected,  and  there 
remains 


c  rC*14>  ay()>  • 

Conformably  we  may  use  the  identification  (129b)  , 
y(1)  =  x,  +0(tr2)  , 

so  that 

f(*>  =  f1  +  0(c“s)  , 

R(14  34)  *  Rt4j4  +  0(C  2)  • 


(140») 
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Thus  (140')  reduces  in  the  linear  approximation  to 


o3r,4m  ■ '  tr  •  (141) 

This  equation  expresses  certain  components  of  the  Riemanniun  curvature 
tensor  in  terms  of  the  measurable  gradient  of  the  3-force  f  .  By  (114),  which 
reads 


_  *2V 

dx,  Bxj  * 

these  components  are  connected  with  the  second  derivative.:  of  the  gravitational 
potential  V.  Thus  we  obtain  by  comparing  these  two  equations: 


a3v  af1 

axj  axj  ax3 


(142) 


This  equation  is  identical  with  (71)  because  of  (104),  but  was  obtained  by 
a  strictly  relativistic  argument.  This  result  is  of  importance  because  it  pro¬ 
vides  the  final  justification  of  the  approach  through  classical  mechanics.  The 
effect  of  the  inertial  forces  is  indeed  removed  if  second  derivatives  of  the 
potential  are  measured  by  means  of  a  gyroscopically  stabilized  apparatus. 


2.3.5.  Relativistic  Conclusions 

The  analysis  of  our  problem  by  means  of  the  general  theory  of  relativity 
has  confirmed  the  result  of  classical  mechanics  that  gravitational  and  inertial 
forces  can  in  fact  be  separated,  at  least  to  a  very  good  approximation,  in  the 
second  and  higher  derivatives  of  the  potential. 
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On  the  one  hand,  this  is  to  be  expected  because  Einstein's  theory  of 
gravitation  goes  over,  as  a  limit  for  a  weak  gravitational  field  and  low  veloc¬ 
ities,  into  Newton's  theory.  On  the  other  hand,  there  seems  to  be  a  disagree¬ 
ment  with  basic  principles  of  general  relativity,  in  particular  with  the  Principle 
of  Equivalence,  according  to  which  gravitational  and  inertial  forces  are  basically 
identical.  The  obvious  conclusion  from  this  principle  is  that  these  two  kinds  of 
forces  are  indistinguishable  and  cannot,  therefore,  be  separated. 

Our  results  are  also  in  apparent  disagreement  with  the  Principle  of 
General  Covariance,  according  to  which  all  coordinate  systems  (inertial  and 
noninertial,  unaccelerated  and  accelerated,  nonrotating  and  rotating)  are 
equivalent  in  principle;  there  are  no  "privileged"  coordinate  systems  (such  as 
inertial  systems)  when  gravitation  is  present.  How  can  we  explain  these  con¬ 
tradictions? 

Both  the  Principle  of  Equivalence  and  the  Principle  of  General  Covariance 
have  played  a  fundamental  heuristic  role  in  Einstein's  considerations  leading  to 
his  theory  of  gravitation  around  1915,  because  these  principles  provide  a  natural 
and  mathematically  obvious  transition  from  the  flat  space-time  of  special 
relativity  to  the  curved  space-time  of  general  relativity.  Since  Einstein's 
heuristic  procedure,  although  it  is  not  a  logical  deduction,  possesses  great 
intuitive  force,  it  is  strongly  emphasized  in  almost  every  introductory  textbook 
on  general  relativity,  subtler  distinctions  being  usually  disregarded  or  bypassed. 


54 


Our  present  problem,  however,  requires  precisely  these  subtler  dis¬ 
tinctions.  We  shall  briefly  present  our  point,  which  is  in  essential  conformity 
with  Fock  (1959)  and  Synge  (1960). 

We  start  with  the  Principle  of  General  Covariance.  The  application 
of  the  relativistic  theory  of  gravitation  to  the  region  of  our  solar  system 
requires  boundary  conditions  at  infinity:  with  increasing  distance  from  the 
attracting  masses  the  effect  of  gravitation  vanishes,  and  the  curved  space- 
time  becomes  flat  at  infinity.  This  fact  permits  the  introduction  of  uniquely 
defined  privileged  coordinate  systems,  the  harmonic  coordinate  systems. 
which  to  a  certain  extent  correspond  to  the  inertial  coordinate  systems  of 
special  relativity  (Fock,  1959). 

Harmonic  (or  "isothermic")  coordinates  date  back  to  de  Donder  in 
1921  and  Lanczos  in  1922;  their  geometrical  and  physical  significance  is 
explained  in  (Darmois,  1927)  and  (Levi-Civita,  1950),  and  their  uniqueness 
in  connection  with  the  boundary  conditions  at  infinity  has  been  investigated  by 
Fock  (1959),  who  used  them  consistently  and  strongly  emphasized  (and  perhaps 
overemphasized)  their  significance. 

General  covariance  is  an  extremely  efficient  means  of  expressing  physical 
laws  as  general  theorems  of  the  geometry  of  space-time  by  means  of  the  tensor 
calculus.  It  does  not  really  assert,  however,  that  in  some  cases  physical  laws 
may  not  assume  a  simpler  form  in  certain  privileged  coordinate  systems.  The 
gravitational  equations  do  assume  a  simpler  form  when  a  harmonic  coordinate 
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system  is  used;  for  instance,  such  a  "quasi- inertial”  system  underlies  the 
expansion  (82). 

As  to  the  Principle  of  Equivalence,  it  is  undoubtedly  true  (because  of 
the  identity  of  gravitational  and  inertial  mass)  as  long  as  the  force  acting  at 
one  point  only  is  considered.  As  soon  as  we  consider  a  region  in  space,  even 
an  arbitrarily  small  one,  however,  we  have  an  objective  criterion  as  to  whether 
a  true  gravitational  field  is  present  or  not.  For  then  we  can  form  the  Riemannian 
tensor  by  means  of  the  forces  and  their  first  derivatives  (the  derivatives  being 
obtained  if  the  forces  are  known  in  an  arbitrarily  small  region);  cf.  equation 
(138).  If  the  Riemannian  tensor  vanishes,  there  is  no  gravitation;  if  it  is 
different  from  zero,  there  is  a  true  gravitational  field.  For  careful  exposi¬ 
tions  of  the  Principle  of  Equivalence  see  (Eddington,  1924,  pp.  39-41)  and 
(Bergmann,  1962,  pp.  204-206). 

The  Riemannian  curvature  provides  a  criterion  for  the  presence  of  a 
gravitational  field,  but  not  yet  a  means  for  the  separation  of  gravitational  and 
inertial  effects,  hi  flat  space-time  inertial  forces  have  an  objective  significance 
because  they  are  due  to  the  deviation  of  the  observer's  coordinate  system  from 
an  inertial  system,  which  represents  a  Cartesian  coordinate  system  in  space- 
time.  hi  the  presence  of  a  gravitational  field,  however,  space-time  is 
irregularly  curved,  and  the  use  of  global  Cartesian  coordinates  is  no  longer 
possible. 
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Thus  a  separation  of  gravitational  and  inertial  effects  is  feasible  only  if 
we  succeed  in  introducing  a  privileged  coordinate  system  sharing  many  desirable 
properties  of  Cartesian  coordinates ,  and  this  separation  is  possible  only  to  the 
extent  to  which  the  privileged  coordinate  system  has  the  properties  of  Cartesian 
coordinates  (it  cannot  have  all  of  them  because  space-time  is  irregularly  curved). 
Probably  the  best  choice  for  such  privileged  "quasi-inertial"  coordinate  systems 
are  the  harmonic  coordinate  systems  mentioned  above.  Inertial  forces,  by 
definition,  will  then  occur  when  the  observer's  coordinate  system  differs  from 
a  harmonic  system. 

This  is  particularly  evident  in  the  linear  approximation,  because  then 
the  gravitational  field  can  be  assumed  to  be  embedded  in  a  flat  space,  and  the 
harmonic  coordinate  systems  can  be  identified  with  the  inertial  coordinate 
systems  of  this  auxiliary  flat  space.  In  fact,  gravitational  and  inertial  force 
fields  are  superposed  linearly  without  interaction;  compare  equations  (77) 

(flat  space),  (80)  (pure  inertial  force  field),  (85)  (pure  gravitational  field), 
and  (86)  (combined  field  of  gravitational  and  inertial  forces). 

Within  the  linear  approximation,  inertial  systems  (or  rather  three 
axes  parallel  to  inertial  axes  in  space)  may  be  defined  even  locally;  they 
coincide  with  global  systems  to  an  accuracy  of  0(cr2),  as  we  have  seen  in 
sec.  2.3.3  ;  according  to  Eddington  (1924,  p.  99),  a  local  terrestrial  inertial 
coordinate  system  will  deviate  from  a  global  inertial  system  by  about  2  seconds 
of  arc  in  a  century. 
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Whereas  geodesics,  in  general,  by  no  means  appear  as  straight  lines  in 
three-dimensional  space  (elliptical  planetary  orbits  are  geodesics  ii:  four-dimen¬ 
sional  space-time!),  it  may  be  shown  that  null  geodesics  representing  light  rays 
in  vacuum  do  appear  as  three-dimensional  straight  lines  (to  0(c~a);  the  bending 
of  a  light  ray  grazing  the  limb  of  the  sun  is  a  higher-order  effect) . 

Also,  by  (93),  "proper  time"  s  coincides  to  0(cfa)  with  "coordinate  time" 
t.  Hence,  even  in  the  "general  theory  of  relativity"  there  are  at  least  three 
"absolute"  concepts  in  the  linear  approximation:  time,  straight  lines  in  space 
(represented  by  light  rays),  and  rotation;  these  concepts  are  absolute  in  the 
sense  that  they  are  practically  the  same  whether  viewed  globally  or  locally,  so 
that  they  can  be  determined  even  by  local  measurements. 

The  deeper  reason  for  the  existence  of  the  linear  approximation  is  that 
gravitational  and  inertial  effects  show  a  characteristic  hierarchy  of  orders  of 
magnitude,  the  differences  being  of  the  order  of  cra  and  therefore  very  marked. 

We  may  summarize  the  three  main  reasons  for  the  possibility  of  separation 
of  gravitational  and  inertial  effects  as  follows : 

1.  The  Riemannian  curvature  tensor  provides  a  means  of  locally 
detecting  the  presence  of  a  gravitational  field. 

2.  Harmonic  coordinate  systems  serve  as  "quasi-inertial"  coordinate 
systems,  with  respect  to  which  inertial  forces  may  be  defined. 

3.  Sharply  distinguished  orders  of  magnitude  greatly  facilitate  the 
separation. 
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The  linear  approximation,  which  is  fully  sufficient  for  our  present  purpose, 
makes  our  separation  practically  unambiguous;  however,  in  more  exotic  gravity 
fields  considered  in  general  relativity,  where  harmonic  coordinates  that  go  over 
into  inertial  Cartesian  coordinates  at  infinity  may  not  even  exists,  a  separation 
of  gravitational  and  inertial  effects  may  well  be  meaningless. 


3.  General  Conclusions 

The  preceding  analysis  has  shown  that  gravitational  effects  can  be 
separated  from  inertial  disturbances  in  the  second  derivatives  of  the  potential 
in  case  of  inertial  stabilization,  and  in  the  third  derivatives  in  the  absence  of 
stabilization.  This  separation  is  possible  because  of  the  locally  different 
behaviour  of  gravitational  and  inertial  fields. 

hi  the  first  derivatives,  forming  the  force  vector,  such  a  separation  by 
means  of  mechanical  principles  is  a  priori  impossible  because  of  the  identity 
of  inertial  and  gravitational  mass.  (It  is,  however,  possible  to  remove  inertial 
disturbances  to  a  large  extent  by  a  suitable  statistical  filtering.) 

Thus,  with  the  advent  of  measuring  techniques  using  equipment  in 
motion,  the  higher  derivatives  of  the  gravitational  potential  may  well  become 
of  increasing  use  to  geodesists.1  Therefore,  we  have  discussed  some  relevant 
geodetic  aspects  in  Chapter  1. 

It  should  be  carefully  kept  in  mind  that  the  gravitational  potential  cannot 
be  determined  without  external  Information  as  to  position,  etc. ,  even  if  the 
second  derivatives  can  be  obtained  directly.  The  reason  is  that  integrations  such 
as  (5)  and  (6)  or  (32)  require  the  knowledge  of  the  position  of  the  measuring  instrument, 
for  instance  of  the  form  (4).  This  position  cannot  be  determined  by  inertial  naviga¬ 
tion  only  because  inertial  navigation  systems  cannot  discriminate  between  gravitational 

1 A  measuring  technique  is  discussed  in  (Thompson,  Bock,  and  Savet,  1965). 
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and  inertial  accelerations,  which  are  the  input;  we  would  thus  be  led  to  a  vicious 
circle  involving  a  bootstrap  operation  as  mentioned  in  the  Introduction. 

hi  the  theoretical  analysis  of  Chapter  2  we  have  restricted  ourselves 
to  the  simplest  case  of  inertial  stabilization,  in  which  the  instrument  axes 
are  kept  parallel  to  axes  fixed  in  inertial  space.  The  results  as  to  separability 
in  the  second  derivatives  hold,  of  course,  also  for  he  case  in  which  the  instru¬ 
mental  axes  are  stabilized  in  such  a  way  as  to  have  a  prescribed  orientation  with 
respect  to  the  earth  (in  this  case,  the  matrix  w,  3  is  different  from  zero  but  given, 
so  that  the  right-hand  side  of  (68)  is  known). 

The  relativistic  analysis  of  the  separability  of  gravitation  and  inertia 
has  provided  an  important  confirmation  of  the  results  obtained  through  classical 
mechanics.  This  analysis  was  necessary  to  remove  the  doubts  raised  by 
Einstein's  Principle  of  Equivalence. 

The  general  theory  of  relativity  is  commonly  accepted  as  the  best  theory 
of  gravitation,  and  it  is  therefore  natural  that  geodesy,  which  uses  gravitation 
to  a  large  extent,  should  be  interested  in  it.  Furthermore,  it  has  a  particular 
appeal  to  geodesists  because  the  geometric  principle  of  Einstein's  theory  of 
gravitation  is  an  extension  of  Gauss'  surface  theory,  familiar  to  every  geodesist, 
to  four -dimensional  space-time,  and  the  relevant  mathematical  technique,  tensor 
calculus,  has  been  introduced  into  geodesy  already  some  time  ago  by  the  work  of 
Marussi  and  Hotine. 
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General  relativity  has  so  far  been  widely  applied  to  problems  of 
astronomy,  and  it  similarly  lends  itself  to  application  to  kinematictil  geodesy; 
as  an  example,  consider  the  "kinematic  gravity  potentials"  (sec.  2.3.1  ) 
which  generalize  the  gravity  potential  in  a  way  necessary  to  take  account  of 
inertial  forces  acting  on  moving  systems.  Matters  may  be  somewhat  com¬ 
parable  to  the  use  of  Hamiltonian  methods  in  satellite  geodesy:  although 
these  methods  are  considerably  more  complicated  mathematically  than  the 
simple  Newtonian  theory,  they  are  much  more  efficient  for  the  solution  of 
some  practical  problems. 

The  foundation  for  an  application  of  the  general  theory  of  relativity  to 
geodetic  problems  may  be  found  in  (Synge,  19 GO);  the  tendency,  the  mathe¬ 
matical  methods,  and  the  applications  are  particularly  appealing  to  astronomers 
and  geodesists.  As  an  example  we  mention  Synge's  "relativistically  valid  geodetic 
survey, "  which  we  have  used  in  sec.  2.3.4. 

Relativistic  methods  should  be  useful  in  problems  such  as  the  kinematical 
measurement  of  gravitational  effects  or  the  interaction  of  inex*tial  positioning 
systems  with  the  gravity  field,  with  which  geodesy  is  likely  to  be  confronted  in 
the  near  future. 
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